A Stochastic Energy Budget Model Using Physically Based Red 

Noise 



M. Weniger, A. Bovier, A. Hense 
January 21, 2013 



Abstract 

A method to describe unresolved processes in meteorological models by physically based 
stochastic processes (SP) is proposed by the example of an energy budget model (EBM). Contrary 
to the common approach using additive white noise, a suitable variable within the model is chosen 
to be represented by a SP. Spectral analysis of ice core time series shows a red noise character of the 
underlying fluctuations. Fitting Ornstein Uhlenbeck processes to the observed spectrum defines the 
parameters for the stochastic dynamic model (SDM). Numerical simulations for different sets of ice 
core data lead to three sets of strongly differing systems. Pathwise, statistical and spectral analysis 
of these models show the importance of carefully choosing suitable stochastic terms in order to get 
a physically meaningful SDM. 



1 Introduction 



The atmosphere itself (or the climate system in 
a wider view) can be seen as a strongly non- 
linear and infinite-dimensional dynamic system 
interacting via a multitude of different time- and 
space-scales of the various flow characterizing 
variables. For numerical simulations a model sys- 
tem of the atmosphere or the climate system can 
treat only a finite number of degrees-of-freedom 
(dof). The effect of the truncated dof on the 
resolved ones can in general not be neglected due 
to the scale interactions. Therefore it is necessary 
to parametrize the effects of the unresolved scales 
on the resolved ones. In principle the stochastic 
character of the unresolved scales has to be 
taken into account. But due to historic reasons 
most parametrization of subscale processes are 
deterministic, they model conditional expectation 
values of moments (mostly second moments) 
given the resolved scales (von Storch 1997). A 
possible problem of numerical weather prediction 
(NWP) and climate modelling using deterministic 
parametrization of subscale/unresolved processes 
is therefore the incomplete consideration of 
interactions between the resolved and sub-grid 
scale processes. A stochastic treatment of these 
parameterizations bears the potential of im- 
proving the simulations and providing a better 
understanding of the stochastic characteristics 
of the simulated atmospheric or climate system 



variable. Furthermore a stochastic model yields 
the natural way to describe the state of a chaotic 
system, i.e. in form of probability densities. This 
in turn is closely connected to the analysis of 
tipping points and regime switches as well as 
occurrences of extreme events e.g. as revealed in 
classical statistical physics (Honerkamp 2002). 
The idea of using stochastic climate models was 
first introduced by Edward Lorenz (1976): "/ 
believe that the ultimate climatic models [...] will 
be stochastic, i.e., random numbers will appear 
somewhere in the time derivatives.^^ However in 
the same article he made a note of caution: "// 
we are truly careful in introducing our random 
numbers, we can likewise assure ourselves that 
the probability of producing an ice age, when 
one ought not to form, is some infinitesimally 
small number''^ . In the same year Hasselmann 
(1976) published a pathbreaking article which 
contains the first mathematical formulation 
of climate models treating weather effects as 
random forcing terms. Since then there have 
been various approaches to incorporate stochastic 
techniques in existing numerical models for 
atmosphere and climate. In many cases this 
was proposed or done without the appropriate 
analyses concerning both the numerical imple- 
mentation and the mathematical and physical 
properties of the stochastic variables or processes. 
While the numerical flaws are easy to avoid in 
self-contained conceptual models when handled 
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with proper care, it is a non-trivial challenge 
to implement a stochastic parametrization in 
an existing more complex model, e.g. a global 
circulation model (GCM), which consists of 
deterministic integration schemes. The second 
aspect is more subtle but equally important: 
the occurrence and the character of a stochas- 
tic parametrization have to be physically justified. 

The importance of carefully choosing the 
stochastic terms is the subject of this article. Ex- 
amplary a energy budget model by Budyko (1969) 
and Sellers (1969) for the global mean temper- 
ature will be considered. Instigated by Hassel- 
manns paper (Hasselmann 1976) it was Fraed- 
erich, who approached this problem by using 
stochastic forcing in form of an additive white 
noise term (Fraedrich 1978). Since then it has 
been an active field of research, see e.g. the review 
papers (North 1981, Ghil 1981), closely connected 
to the phenomenon of stochastic resonance (Benzi 
et al. 1982). 

After defining the deterministic framework (|2.ip 
we will point out problems of the white noise 
ansatz regarding the model physics (|2.2.ip . A dif- 
ferent, approach using physically based stochas- 
tics is proposed with a focus on spectral ice core 
data analysis (|2.2.2p . It turns out that different 
data sets lead to three parameter classes for the 
stochastic terms. For each class one examplary 
model is considered. The simple mathematical 
structure of the model exhibits very effective nu- 
merical simulations for nonlinear stochastic terms 
(j2.2.3|) which in turn allows comprehensive statis- 
tical evaluation (jSj. Particular emphasis is placed 
on the pathwise analysis of the models, since basic 
statistical properties like mean value and variance 
are unable to capture the crucial characteristics, 
i.e. bifurcation structure, timescales of crossings 
between stable states and correlation between in- 
solation and temperature. After getting an idea of 
the models behaviour by looking at single realiza- 
tions p.ip . marginal distributions of the temper- 
ature at local insolation extrema (|3.2p as well as 
coherence analysis in the frequency domain p.3p 
yield solid evidence on the differences between 
these models. Depending on the choise of ice core 
data we get fundamentally differing temperature 
processes. Since there is no reason why a par- 
ticular data set should be preferred, we come to 
the conclusion that the model is inapt to be fitted 
on the type of data available, which is discussed 
in detail in p.4p . This however emphasizes the 
importance of choosing the stochastic terms very 



careful even in such a simple case. The main focus 
of this study are not the results of the stochastic 
EBM but rather the method used to derive them, 
i.e. 

1. Identify a variable suited to be describes 
stochastically 

2. Use spectral analysis of data to determine the 
distribution of the SP 

3. A rigorous numerical implementation, which 
will be a nontrivial task in more complex 
models involving stochastic partial differen- 
tial equations 

4. Statistical and pathwise analysis of the 
stochastic dynamic model 

1.1 Terminology 

Looking at various publications regarding stochas- 
tic models there is no consistent terminology. 
We will understand "Stochastics" as the field of 
mathematics which summarizes "Statistics" and 
"Probability Theory". Where the later includes, 
the branches "Stochastic Processes" and "Ran- 
dom Variables". The term "Stochastic Parame- 
terization" will be used to describe models where 
state independent parameters are replaced by ran- 
dom variables. Models in which time- or state- 
depending stochastic processes are used will be 
denoted as "Stochastic Dynamic Models" (SDM). 
On the one hand this emphasizes the fact that the 
stochastic term itself is part of the systems dy- 
namic and on the other hand it draws the connec- 
tion to the mathematical field "Stochastic Dynam- 
ics" . Linked to a distinct terminology but with 
direct mathematical implications is the differen- 
tiation of the occuring types of differential equa- 
tions. Following (Jentzen and Kloeden 2009b) we 
will distinguish between the following types of dif- 
ferential equations 

• deterministic: ordinary differential equations 
(ODE) / partial differential equations (PDE) 

• random: RODE / RPDE 

• stochastic: SODE / SPDE 

Random differential equations (RODE/RPDE) 
are ODEs/PDEs with a stochastic process Wt 
in their vector field, whereas stochastic differen- 
tial equations (SODE/SPDE) can be character- 
ized by the appereance of differential terms dWt 
of a stochastic process, i.e. the unknown quan- 
tity itself is a stochastic process. Especially these 
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two classes are often mixed up in meteorological 
literature. They show essential mathematical dif- 
ferences: RODE/RPDE can be interpreted and 
analyzed pathwise as their deterministic counter- 
parts. Note that random variables, which can be 
seen as time constant SPs, are included in the 
definition of RODE/RPDE. For time dependent 
SPs however the emerging vector field is usually 
not differentiable because most driving SPs are 
only Hocldcr continuous in time. This doesn't 
prohibit the use of the deterministic calculus but 
leads to lower convergence orders for the numeri- 
cal schemes (Kloedcn and Jentzcn 2007). 
SODE/SPDE require a different kind of calculus. 
The vast majority of mathematical literature re- 
garding stochastic differential equations is based 
on the Ito calculus (Ito 1951), which requires a 
numerical implementation different from the de- 
terministic one. While techniques for SODEs are 
well developed, e.g. (Kloeden and Platen 1992), 
"t/ie numerical solution of SPDEs is at a stage 
of development roughly similar to that of SODEs 
in the 1970s'' (Jentzen and Kloeden 2009a). In 
this article we will have to deal with RODEs and 
SODEs. 



ical processes which may have an effect on the 
model variables. For this we take a look at each 
occuring physical variable whose dynamics are not 
captured by the model, especially the ones which 
are approximated by a constant value. For two 
kinds of physical values a stochastic representa- 
tion is uncalled-for: "true" constants, e.g. Stefan- 
Boltzmann-constant, and non constant variables 
which could accuratly be described in a determin- 
istic way, e.g. incoming solar radiation depending 
on Milankovich cycles. The interesting ones are in- 
exact approximations of processes which cannot be 
described in a deterministic way within the models 
framework. Of course stochastic representations 
arc also limited by the models space- and time- 
scales, which provides another criterion in the se- 
lection process. When a eligible variable (or a set 
of eligible variables) is found we need some kind 
of input on the distribution of the stochastic proc- 
cess. This may be information based on models, 
theoretical insights or data. Deducing the com- 
plete distribution of the SP can be a nontrivial 
task and should be done carefully since even small 
variations may lead to significant changes in the 
dynamic structure of the system. 



1.2 A Method for Physically Based 
Stochastic Models 

In most cases the starting point for the develop- 
ment of a stochastic model is a deterministic one 
which captures the gouverning dynamics. The role 
of stochastic processes is to describe the dynamics 
of fluctuations and/or uncertainty of model vari- 
ables and/or sub grid processes with an emphasis 
on the word "dynamics": whereas stochastic pa- 
rameterizations can be used to represent uncer- 
tainties of model parameters, initial values and 
other time independend variables, SP provide a 
tool to lay a hand on time and state dependend in- 
teractions between subgrid fluctuations and model 
variables. To refer to the early meaning of stochas- 
tic as the "art of guessing" it is indispensable to 
carefully choose the correct SP for a given model. 
In a more recent article Penland and Ewald (2008) 
stated "Simply replacing the fast term with a 
Gaussian random deviate with standard deviation 
equal to that of the variable to be approximated, 
and then using deterministic numerical integra- 
tion schemes, is a recipe for disaster". While the 
use of a correct stochastic numerical approach is 
equally important this paper primarily deals with 
the problem of finding a physically based SP. As 
the first step we have to identify unresolved phys- 



2 An Energy Budget Model 

2.1 The Deterministic Framework 

The deterministic EBM is given by an ODE for the 
global mean temperature Tt characterizing the ra- 
diation balance between net short wave radiation 
and outgoing long wave black body radiation 

c dTt = [St{l- a{T)) - a{Tt ~ ^Tf]At. 

Here c = 3 * l(f'J/{m'^K) denotes the global 
heat capacity, S the insolation, a{T) the albedo, 
cr = 5.67* 10~^iy/(m2ii:4) the Stefan-Boltzmann- 
constant and AT = 32.6i^ the temperature offset 
between surface and the top of the atmosphere. 
Since we are dealing with a model with spatial 
dimension zero, obliquity and precession are not 
accounted for and the solar forcing is solely depen- 
dent on the eccentricity of the earth orbit (FigH]). 
Note that the absolute value of S only fiuctuates 
in the small interval between ZAl.hK and 342. Oif. 
The albedo was parameterized by deriving the 
values needed to reproduce stable equilibria at 
typical ice- and warm-age temperatures, which 
were estimated from ice core data, and taking into 
account the present state: (Ti = 278.9/!r, a(ri) = 
0.3901), (T2 = 288.0A',a(r2) = 0.2950), (T3 = 
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290.3if, a(r3) = 0.2676). Between these points a 
stepwise linear interpolation was used (FigH]). 
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Figure 1: Albedo parameterization and insola- 
tion cycle used in the deterministic 
framework 



Looking at the dynamic structure of this system 
we observe two stable equilibria at Tice = 278.9i4r 
and Twarm = 290.3K, which are very unsuscep- 
tible to variations in insolation. The tempera- 
ture value of the third equilibrium strongly de- 
pends on the solar forcing and can be found at 
T e (283.57^:, 287.3if ). Note that if the insolation 
would vary in a slightly broader interval, we could 
observe areas where this equilibrium wouldn't ex- 
ist, which is why it will become a meta-stable state 
in the stochastic model. For an arbitrary initial 
temperature value the system will eventually set- 
tle in one of the stable states without the possi- 
bilty to cross between these states due to varying 
insolation. 



2.2 A Stochastic EBM 

2.2.1 Motivation 

Looking at ice core data immediately reveals, that 
the temperature is a strongly fluctuating process 
varying between ice- and warm-ages. The essen- 
tial missing trait of the deterministic model is 
the ability to switch between its steady states. 
One common approach is the use of stochastic 
forcing in form of an additive white noise term 
(Fraedrich 1978). While being a reasonable first 
step this ansatz is disputable concerning 

• additivity 

There is no temperature independent vari- 
able in the model which leads to the question 
what physical process should be described by 
a temperature independent noise term. 

• whiteness 

Since the physical processes described by this 
model acts on large space- and timescales it 
is hard to imagine why a non-correlated noise 
should adequately capture the occurring fluc- 
tuations. 

This was the motivation to use a different ap- 
proach. Rather than using a purely mathemat- 
ical tool to modify the model in order to pro- 
duce the desired effect, we looked at the variables 
which are already used by the deterministic EBM 
to describe one of them as a physically based SP. 
For a successful stochastic formulation a variable 
has to satisfy the following conditions: foremost 
there have to be significant fiuctuations of a global 
mean because the model doesn't resolve spatial 
dimensions. These fluctuations have to act on a 
timescale which can be covered by the model, i.e. 
not smaller than one year. And finally we need 
evidence on the distribution of the SP. This can 
be achieved by means of another model resolving 
the variable in question or ideally through data 
analysis. 

2.2.2 Data Analysis 

Considering these demands we chose AT and 
its variations due to non constant concentra- 
tions of greenhouse gas. We used the follow- 
ing sets of ice core data: Vostok (Wahlen et al. 
2000), DomeFuji (Kawamura et al. 2007), Byrd 
(Ahn and Brook 2008), EPICADome (Monnin 
et al. 2004), TaylorDome(Indermuehle et al. 1999), 
TaylorDome2 (Indermuehle et al. 2000). A graph 
of the Byrd data (Figl3|) shows examplary that the 
scaledemands are met and provides direct input 
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Figure 2: Byrd ice core data 



regarding the distribution of the SP. One has to 
be careful on choosing the right data here because 
there is a strong correlation between CO2 concen- 
tration and temperature which involves the danger 
of predefining the models equilibrium structure. 
In this simple model we would like to have a SP 
which describes the "natural" small fluctuations of 
AT excluding feedback effects, i.e. a SP indepen- 
dent of the model variable T, and analyze its ef- 
fects on the systems dynamic structure. Thus the 
perfect data would be a large sample of CO2 data 
on a time interval where the temperature stays 
constant. Unfortunately this kind of data is not 
available but we can come close by analyzing de- 
trended dimensionless time-series without cross- 
ings between ice- and warm-ages. The first step is 
the conversion of CO2 data to temperature anoma- 
lies, where a doubling in CO2 concentration yields 
a AT increase of 2K and we get 



AT' 



log 2 



log 1 



dctrcnd(C02) 



AT' is the desired data of temperature fluctu- 
ations, detrend(C02) the lineary detrended data 
from ice core drillings and (CO2) the mean CO2 
level. In order to characterize the SP wc did a 
spectral analyzation of AT' time-series which in- 
dicates a red noise distribution (Fig 13]). In math- 
ematics this AR(1) process is known as Ornstein- 
Uhlenbeck process (OUP), or Gauss Markov pro- 
cess, and can be defined by the Ito SODE 

dATt = -e{ATt - fi)dt + \/Z) dWt. 

In general the decision between Ito and 
Stratonovich calculus is an important one, 
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Figure 3: Spectral analysis of Byrd, Taylor and 
Vostok ice core data 



in this special case however it doesn't make a 
difference because the noise term is not state de- 
pendent. The OUP is defined by three parameters 
fi, Q and D. Since we used dimensionless time- 
series our approach uses a constant expectation 
value fi which is the AT used in the deterministic 
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model. The remaining parameters, which are 
responsible for the fluctuations, were fitted to 
the spectral data and confirm the validity of the 
red noise assumption. Note however that the 
confidence interval for the parameter values of 
each time-series is quite large due to low data 
size. In addition the parameters calculated from 
different time-series differ strongly and it is not 
obvious which parameters we should trust: 
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These time-series can be divided into three 
classes [DomeFuji and Vostok], [Byrd, EPI- 
CADomc and TaylorDome2] and [TaylorDome]. 
There are pro and contra arguments for each 
class: the first one consists of very large data 
sets covering long time periods. While the high 
amount of data points is certainly a positive as- 
pect it implies the problem of embedded crossings 
between ice- and warm-ages. The second class 
contains the majority of time-series where each 
member yields similar parameters. The third 
class consists of only one scarce time-serie with 69 
data points. Nevertheless it leads to the weakest 
noise signal, so one would be "on the safe side" 
using these parameters. Since there is no good 
argument to exclude either class we analyzed the 
system for one time-series of each category: Byrd, 
TaylorDome and Vostok. 

2.2.3 Numerical Aspects 

The stochastic EBM is a twodimensional system 
containing one RODE for Tt and one SODE for 

c dTt = [Stil - a{Tt)) - <j{Tt - ATt)^]dt 

dATt = -e(ATt ~ ^i)dt + VD dWt. 

The second equation can be solved analytically 
yielding 

ATt = AToe-®* + //(I - e"®*) 




This allows the process ATt to be simulated in 
a computational cost effective way by generating 



Gaussian random variables with a variance corre- 
sponding to the exponential time scale e^^*. The 
remaining RODE can be solved by common deter- 
ministic numerical schemes, allthough they con- 
verge at a slower rate, since the process ATt is not 
differentiable but only Hoelder continous (Klocden 
and Jentzen 2007). We used a fourth order Runge- 
Kutta scheme with a timestep of one year. 

3 Results 

3.1 Samplepaths 

Looking at a samplepaths (Fig|4]) of the three sys- 
tems gives a first intuitive impression of the models 
characteristics: 

• there are at least two stable equilibria coin- 
ciding with the deterministic ones 

• the system can perform jumps between these 
equilibria 

• the models differ strongly regarding decorre- 
lationstimes of Tt and strength of correlation 
between Tt and St 

The samplepaths suggests a correlation between 
solar forcing and temperature depending on the 
strength of the SP. The coupling of solar and 
stochastic forcing with the temperature can be 
explained on a heuristic level, i.e. ignoring the 
systems inertia, by looking at schematic potential 
plots (FiglS)) . Examplary we will consider positive 
fiuctuations of ATt where the system (illustrated 
by a black ball) is in an ice-age. Positive fluctua- 
tions of insolation St and offset temperature ATt 
cause the potential to tilt towards the warm sta- 
ble state. Since the insolation changes alone are 
not capable of causing crossings between ice- and 
warm-ages we have to take a closer look at AT. 
Qualitativly the values of AT can be broken down 
into three categories: 

• an area around the mean value 

• an one-sided area further away from the 
mean. i.e. a "weak" extremum 

• an one-sided area far away from the mean, i.e. 
a "strong" extremum 

While AT occupies the first area (light gray) 
the system stays in its current state unsusceptible 
to the state of solar forcing, i.e. behaves just 
like the deterministic system. When AT adopts 
a weak extremum (medium gray) the system 
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Figure 4: Samplepaths for temperature (black) and insolation (grey) of three models based on different 
ice core data. 
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follows the course of insolation, e.g. it jumps 
from an ice- to a warm-age if and only if the 
insolation is strong at that point in time. In this 
case the stochastic forcing amplifies the solar 
forcing but has no direct effect on the state of the 
system. This changes when AT reaches a strong 
extremum (dark gray): the system will jump 
regardless of the insolation level, e.g. it will jump 
from an ice- to a warm- age even if the insolation 
is minimal. The stochastic forcing is strong 
enough to cause crossings between stable states 
even against the effect of the solar forcing. Since 
the OUP is a Gaussian process, the probabilty 
for extremes corresponds to e"^^"^'"^^"^'^-' . For 
a system with small fluctuations in AT(, e.g. 
Taylor with var(ATt) « IQ-^K"^, jumps of the 
first kind, which are following the insolation, are 
much more likely than jumps of the second kind. 
For a system where AT varies strongly, e.g. Byrd 
or Vostok var(ATt) sa 10~^K^^ crossings of the 
second kind become more common. Note that 
the decorrelationtime of T* corresponds to the 
solar forcing in the first case but to the stochastic 
forcing in the second case. This explains the 
different decorrelationtimes comparing the Byrd 
and Vostok model. 



3.2 Marginal Distributions at Local 
Extrema of the Solar Forcing 

Following the intuition gained by samplepaths 
and heuristic potential considerations we take a 
first step to describe statistical characteristics of 
the models by analyzing marginal distributions 
of Tt at local extrema of the solar forcing St- 
We expect the Taylor model to follow the course 
of insolation and show one-sided distributions, 
whereas the Byrd and Vostok models should 
exhibit close to uniform distributions. In Figl6] 
the temperature was divided into three intervals 
corresponding to ice- and warmages and to a 
state inbetween these two, which is occupied by 
the system during a crossing. Examplary four 
local extrema were chosen and the temperature 
was evaluated with a time lag of Ska for 4000 
realizations of each model. This value will be 
derived in the following section 13.31 The Taylor 
model follows our expectations espacially during 
strong insolation extrema. Even in the case of 
weaker local extrema the distributions strongly 
favor the corresponding temperature state. Con- 
trary to the differing samplepaths for the Byrd 
and Vostok systems, the marginal distributions 



show a strong resamblance. On the one hand this 
shows how important a pathwise analysis of a 
stochastic model is, to capture its characteristics, 
on the other hand it is an indication that both 
models have similar correlation between insolation 
and temperature, independend of the differing 
decorrelation times. The later show up in form 
of the differing amount of transition states, which 
hints at more occuring crossings in the Byrd 
model. 



3.3 Coherence 

Of course neither a single path nor marginal dis- 
tributions at fixed points in time are sufficient to 
characterize a stochastic system. Therefore we an- 
alyzed the crosscorrelation and coherence between 
insolation and temperature based on the 4000 re- 
alizations of each model we used in the last step. 
The impression gained thus far are verified by the 
crosscorrelation (Fig I?]), which is much higher in 
the Taylor system than it is in Byrd in Vostok 
systems. A maximal crosscorrelation of about 0.8 
at a time lag of 8ka and the peaks at time lags 
of about 4Q0ka, which relates to the eccentric cy- 
cle, suggests that the Taylor system is dominated 
by its insolation. The same periodic structure can 
be observed for the Byrd and Vostok model but 
with a much smaller amplitude. In these models 
a maximal crosscorrelation of about 0.2 indicates 
only slight influences of the insolation on the tem- 
perature but a dominance of the stochastic forcing. 
The periodic crosscorrelation structure suggests 
to take a look at the frequency domain. The 
(quadratic) spectral coherence of two signals x, y 
is defined by C^^ = \Gxy\'^ /G^xGyy, where 
Gxyif) = {X*{f)Y{f)) is the cross-spectral den- 
sity using the Fourier transformed signals X, Y. 
For ergodic linear systems the coherence is a mea- 
sure for the causality between the input x and the 
output y. Allthough we are not dealing with a 
linear system, the fact that the insolation St is a 
deterministic signal leads to a closely related in- 
terpretation. Denoting the Fourier transforms of 
insolation and temperature as S and T yields 

2 \{s*v\' _ m{T)\' 

'^^ (|5p)(|rp) |5p(|rp) 
(irp) [\{r)\'^') ' 

The coherence between S and T does not depend 
on the insolation S. It takes on values between 
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and 1 and can be interpreted as a measure for the 
deterministicity of the temperature T: a system 
with no fluctuations has var(T) = and therefore 
C|n = 1. If the normahzed variance is large Ct 
goes to zero. The frequency interpretation shows 
the impact of the stochastic forcing on the individ- 
ual timescales. All models show peaks at periodic 
times of lOOfca and 120fca (Figl?]), which coincides 
with the spectrum of the eccentric cycle, i.e. the 
solar forcing. As expected the peaks for the Taylor 
model are close to one, whereas Byrd and Vostok 
models produce significant lower values. 

3.4 Conclusions 

Although the SP were carefully constructed in 
each model, the emerging systems differ funda- 
mcntcly. While the Taylor systems yields the 
most realistic results, there is no solid ground 
to utterly dismiss the other data sets. The 
underlying problem is the discrepance of spatial 
dimensions between data sets and model space. 
Ideally we would need data which allows the 
derivation of global mean fluctuations of AT. 
Another factor are feedback effects between CO2 
levels and temperature which are included in the 
data but not in the model. A third obstacle is 
the large amount of data points one would need 
to determine the SP with high accuracy and 
confidence. When these requirements are met the 
spectral fitting of SPs onto the data provides a 
powerful tool, to include unresolved physics into 
the model. However the diversity of the three 
systems emphasizes the importance of choosing 
the stochastic terms very careful even when 
extending such a simple deterministic model, 
which puts additive white noise approaches into 
question. 
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Figure 6: Marginal distributions at local extrema of the solar forcing St 
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Figure 7: Crosscorrelation and coherence of various ice core data 
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